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Abstract 

We study the non-equilibrium properties of a dynamical fluid composed of 
quasi-particles whose mass is temperature and charge chemical potential 
dependent, in kinetic theory under the relaxation time approximation. In 
particular we calculate the scaling behaviour of bulk viscosity 'C' near the 
QCD chiral phase transition in the 3d 0(2) universality class. It is found 
that the bulk viscosity 'C does not show a divergent behaviour near the 
QCD chiral phase transition. This scaling behaviour of 'C' prevails in the 
presence of Gold-stone modes that arise due to the explicit breaking of 
continuous 0(4) symmetry. On contrary these modes have a significant 
effect on the scaling behaviour of specific heat Cy which diverges at the 
critical temperature Tc- 
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Introduction 



Transport properties of a strongly interacting matter, encoded in coefficients such 
as shear and bulk viscosities, which describe the hydrodynamical response of the 
system to energy and momentum density fluctuations are of particular importance 
for understanding the nature of QCD matter. Their values and properties not only 
carry information on how far the system appears from an ideal hydrodynamics 
but can also provide a relevant insight into the fluid dynamics and its critical 
phenomena [H El E]- Even though tremendous efforts have been put forward from 
both theoretical and experimental sides, their exact determination still remains a 
challenging task, especially in the regime of experimentally accessible temperature 
and densities. 

Recently in an analysis of the non-equilibrium properties of a quasi-particle medium 
in the kinetic theory, under the relaxation time approximation, the scaling be- 
haviour of bulk viscosity ( near the QCD chiral phase transition was studied 
Using the static critical exponents for 3d 0(4) symmetric spin models, it was 
pointed out that the bulk viscosity does not show the singular behaviour near the 
QCD chiral phase transition, in contradiction to the results obtained using the 
Staggered and Wilson fermion formulation of lattice QCD [5][n], wherein a sharp 
rise of C at T = Tc was reported. It was assumed that an extension of relaxation 
time approximation that allows for the proper incorporation of long range fluctu- 
ations of soft modes might result in the divergences of ( at the second order 0(4) 
transition. 

In general the QCD phase transition in the limit of Uf massless quark flavors is 
controlled by SUiinf) x SU[i{nf) chiral symmetry. Renormalization group argu- 
ments suggest that QCD with three degenerate light quark flavors has a first order 
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phase transition [7], whereas the two flavor theory is expected to have a second 
order phase transition. In the latter case the SUl{2) x SUr{2) chiral symmetry 
is isomorphic to 0(4) and the transition therefore is expected to belong to the 
same universality class as three dimensional 0(4) symmetric spin models [8] [9]. 
Depending on the value of strange quark mass the QCD phase transition in the 
limit of vanishing light quark masses may be first order or continuous transition 
still belonging to three-dimensional, 0(4) universality class [7]. 
However none of the above mentioned lattice discretization schemes (Wilson/ Stag- 
gered) for the fermion sector of QCD preserve the full chiral symmetry of QCD 
Lagrangian. At non-zero lattice spacing the Wilson fermion formulation does not 
preserve any continuous symmetry related to the chiral sector of the QCD, while as 
the Staggered fermion formulation preserves at least an 0(2) symmetry, that gives 
rise to the massless Gold-stone modes in the chiral limit. Due to the existence of 
these modes singularities are expected on the whole coexistence line, for the 0{N) 
models with N > 1 and dimensions 2 < d < 4 [TU] [TT], in addition to the known 
singular behaviour at T^. For N=2 and d=3 there exists a rigorous proof [T2] [T3]. 
Keeping into account the above mentioned facts, 1. presence of only 0(2) symme- 
try in case of staggered fermion formulation of lattice QCD, 2. the generation of 
Gold-stone modes due to spontaneous breaking of continuous 0(4) symmetry to 
0(2) and 3. the resulting additional singularities due to these modes, we therefore 
calculate the scaling behaviour of bulk viscosity ( for a quasi-particle medium in 
the 3d 0(2) universality class and evaluate the possible role of the massless Gold- 
stone modes on the scaling behaviour of ( near the QCD chiral phase transition. 

I. Shear and bulk viscosities in a quasi-particle model. 

Let us consider a dynamical fluid composed of quasi-particles where the mass of 
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each quasi-particle is temperature and charge chemical potential dependent. The 
calculation of the transport coefficients bulk viscosity ( and the coefficient of vis- 
cosity T] requires the knowledge of molecular distribution function f{p,x) in terms 
of which a complete description of the macroscopic state of a system is assumed 
to be possible. Under the relaxation time approximation the distribution function 
f{p,x) satisfies the linear differential equation [H] 



ot or m av tq tq 

This equation characterizes the deviation of a system from the state of equilibrium, 
quantified by the equilibrium distribution function /q. Here m and v are the 
mass and velocity of constituents of the medium, that is under the effect of an 
external force F. Here tq is the mean time between successive collisions among the 
constituents of the medium and it is by these collisions a system is assumed to 
restore its equilibrium state. 

For a system that is slightly removed from equilibrium, one can approximate the 
non-equilibrium distribution function f{p,x), by the perturbative expression 



f{p,x) = fo + fi (2) 

where, fi« fo- Under this approximation the left side of Eq.l can be approxi- 
mated by neglecting terms in fi which results in an expression 



^/ = ^ + v ■ — + — ■ — = 3 

at or m av tq 



Here we intend to study a system which is not under the influence of any external 
field, which for-example in our case could be an externally applied magnetic field. 
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therefore we can use, F = in Eq.(3) and we get 



(4) 



Here, = {E, p) is tlie momentum four vector and = {d^, &'') is the four 
operator. Also we have used the relation £J/p =l/v between the energy E — 
•\/p^ + w? and momentum p of a relativistic particle which in absence of external 
force F is having a pure thermal motion. One can use the above calculated variation 
5/, to calculate the variation of energy- momentum tensor from its equilibrium 
value This variation of distribution function 5f can be directly used to calculate 
the variation of energy-momentum tensor from its equilibrium value, that in 

itself contains the information regarding the transport coefficients of the medium. 
For a system composed of fermions and/or bosons that is slightly removed from 
equilibrium, the variation ST^'' can be written as 



where dF = g((i^p/(27r)^ is the integration measure in the momentum space with 
the degeneracy factor g associated with the particle quantum numbers. If we use 
the value of 5f from Eq.(4) in Eq.(5) one has 




(5) 





Here /o and /o are the equilibrium distribution functions for particles and anti- 
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particles respectively and are given by 



/o(/o 



(7) 



where u is the flow velocity and fi is the chemical potential related to any conserved 
charges. The ±1 correspond to fermion and boson statistics, whereas ±/i to particle 
and anti-particle contributions, respectively. If one assumes that the particles and 
anti-particles instead of carrying a fixed value of the mass M, have a dynamical 
mass given by M = M(T, fi), than the above variation ST'^'^ can be further written 
as II] 



+ J rfr^/o(l±/o) g/(p,T,^) 



where 



g;j(p,T,/i) 



p2 fdP\ f dE dE 

2E 



dP\ fdE_ ^ 
dn J \dn 



dkUi6 



kl 



and 



Wki = dkU^ + diu^ - -SkAu' 



using in Eq.(8) the energy conservation / dTE [tq/o (1 ± /o) Qf + to/o (l ± /o) qj] 
and comparing the resulting expression with 



6r^ = -cSijdkU^ - vWij 



(9) 
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where ( and rj are the bulk and shear viscosity of the medium, one obtains follow- 
ing expressions for the transport coefficients of the medium [1] 



1 



IST J 



[grUl ± /o) + -gfMl ± /c 



(10) 



and 



(27r)3~E 



[(7ro/o(l±/o)+^ro/o(l±/o)] 



1 

- - fdP 
— [grUl ± /o) - ^fo/o(l ± /o)] [j^ 



The derivatives dP/de\n and dP/dn\^ which appear in Eq.(ll) can be written in 
terms of various susceptibilities Xxy = d^P/dxdy as follows 



dP\ ^ -sXmm ~ ^X^lT 



(12) 



and 



dP\ 
dn J 



(13) 



Here n = dP/d^ is the net particle number density and s = dP/dT is the 
entropy density. The specifiv heat Cy is given by 



Cv=T 



Xtt — 



XfiT 

XflfM 



(14) 



using the expressions of Eq.(12)- Eq.(14) we will now study the scaling behaviour 
of transport coefficients near the quark- hadron phase transition region. 
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II. Transport coefficients near the QCD chiral phase tran- 
sition. 

To calculate the transport coefficients for a thermodynamic system, one requires 
to work with a specific model for particle interactions. The effect of these in- 
teractions is to change some of the observables of the system for example, the 
particle mass M and charge chemical potential [IS] [IS]! which serve as an input 
in the determination of equilibrium distribution functions /o, /o, of particles and 
anti-particles respectively. These distribution functions together with the ther- 
modynamic variables, for example temperature T and energy E can be used to 
determine the transport coefficients ( and 77 for a given medium. However there 
are some generic properties of bulk viscosity ( that can be studied in a model in- 
dependent way, through the universality arguments. In this context, of particular 
interest is behaviour of ( near the phase transition. 



III. ScaUng analysis of bulk viscosity (. 



To study the scaling of bulk viscosity ( we proceed as follows. First we calculate 
the scaling behaviour of ( in the 3d 0(2) universality class without incorporating 
the symmetry breaking effects that are applicable in view of breakdown of contin- 
uous 0(4) symmetry to 0(2). Then we will incorporate the effect of Gold-stone 
modes in our calculation and evaluate the effect of these modes on the scaling 
behaviour of bulk viscosity near the quark-hadron phase transition region. 
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IV. Scaling of the bulk viscosity ( in the 3d 0(2) uni- 
versality class. 



In the vicinity of the chiral phase transition, the thermodynamic Free energy re- 
ceives contribution from the analytic or regular Gr and non-analytic or singular 
part Gs [T7]. This non-analytic part Gs of free energy results in the singular 
behaviour of different observables near the phase transition [IB] [IS] [20] [IZ]- To 
parametrize the scaling behaviour of Gg near the chiral phase transition let us con- 
sider a general thermodynamic scaling formula for the singular part of free energy 



Here /3 and S are the critical exponents for a given universality class. M is the order 
parameter of the transition and 't' is the reduced temperature variable defined 
as t = t + A /i^ [17J where , i = (T — T^/Tc and Ji = fi/Tc- In the discussion 
to follow we will calculate all the relevant quantities required for examining the 
scaling behaviour of bulk viscosity (, using the scaling form of Gs as mentioned 
above. We will limit our discussion to the /x = case only. 

Using the relation H = dGs/dM , where magnetic field H is a variable conjugate 
to M, together with scaling law for Gg as given by Eq.(15) we obtain 



(15) 



(16) 



Here /' is partial derivative of / with respect to M. Now calculating the inverse of 



9 



Eq.(16) we get the scaling law for M 

M = t^c (Ht-^^) (17) 

where c = f~^. Now we will use the scahng function Gg as defined in Eq.(15) 
together with generallised reduced temperature t = i + Afj^ and calculate all the 
relevant quantities required for examining the scaling behaviour of bulk viscosity 
C,. In the following we will consider only the transition point at // = 0. As already 
mentioned we are dealing with a system which is not under the influence of any 
external force therefore we will use H=0. 

For T — > Tc at zero field H, the scahng behaviour of order parameter M is 

M~t^c(0)~t'^ (19) 



The scaling behaviour of dM/dT is 



dM 



(20) 



Using the scahng behaviour of M and dM / dT, the scaling of susceptibility Xij.ij. — 

d'^Gg/dji'^ can be readily computed. 

Taking double derivative of Gg with respect to n we get 



Now since t = t + Afi'^ = (T - Tc)/Tc + A {fi/Tcf, we have dt/dfx = 2Afi/Tc 
and dH/djj?' — 2A/Tc. Using these expressions together with the scaling of M as 
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in Eq. (19) the scaling form of x^^^ becomes 



l-a 



(21) 



where we have used the relation /3{1 + S) — 2 — a. Next we calculate the scahn^ 
form of xtt — d^Gs/dT^. Using the scaling function Gg as in Eq.(15) and taking 
double derivative with respect to chemical potential /i we get. 



Xtt ^ t^^^^+^^-^^ ^f{Mt-^) + 2f{Mt-^)t-^ + t-^^f"{Mt-^)^{dt/dTf 

+ ^[/3(i+5)-i] r/(Mt-^) + r^/(Mr^)] (dH/dT^) 
+ 1^^^-^^ [/'(Mt-^)l {dM/dT) {dt/dT) 



Using dt/d/i — 2Aii/Tc, together with the scaling law for order parameter M as in 
Eq.(19) we get 



Xtt ^ t' 



(22) 



Now we calculate the scaling of susceptibihty x^t — d^Gs/djidT. Using again the 
scahng function of Gg as in Eq.(15) we get. 



X^,T ^ t^^^^+^^-^^ \f{Mt-^) + 2f{Mt-^)t-^ + f"{Mt-^)t-^P] {dt/dn) {dt/dT) 



^ ^[/9(l+5)-l] 



After substituting the values of different derivatives we get 



(23) 



Therefore the scaling behaviour of specific heat Cy given by Eq.(14) near the phase 
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transition at n=0 becomes 



Cv = T 



Xtt — 



XflfM 



(24) 



Using the scaling laws for x^^, XfiT, Xtt and specific heat Cy, the scaling behaviour 
of bulk viscosity ( is readily calculable. From Eq.(ll) the divergent behaviour of 
( can be written as |1] 



/•div 



M2 



de J dT ^ dfi J \dn J ^ dfi \dn 



(25) 



Using Eq.(12), Eq.(13) together with the scaling laws of X/i/^, Xf^T, Xtt and Cy we 
have 



/ dP\ ^ sx^,^, - nx^,T ^ ^„ 



f dP\ ^ riTxTT + {nil - sT)x^t - sfiXfifi ^ q 



Also we have the scaling relations, dE/dT = {M/E){dM/dT) ~ t^/^-i ^^^^ 
dE/dfi = {M/E){dM /dfi) ~ Therefore the scaling of bulk viscosity ( 

is given by the relation 



(26) 



Now using the critical exponents for 3d 0(2) model, a= -0.01722 and (3= 0.349, 
the singular part of bulk viscosity, C^^'^, vanishes near the phase transition region. 
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Table 1: Critical exponents a, /3, 7, 6 and the universal constant £2 for 
three-dimensional 0(4) and 0(2) universality class \TI\ |23] . The critical 
exponents are connected by the relations 7 = /3(5 — 1) and (3 (1 + 5) = 2 — a. 



Model 


a 


(3 


7 


S 


C2 


0(4) 
0(2) 


-0.21 
-0.01722 


0.380 
0.349 


1.453 
1.319 


4.824 
4.780 


0.666(6) 
0.592(10) 



V. Effect of Gold-stone modes on the scaling of bulk vis- 
cosity (. 

The spontaneous breaking of a continuous 0(N) symmetry results in the generation 
of Gold-stone modes. For an 0{N) invariant nonlinear a- models these modes 
lead to the singular behaviour of susceptibilities near the phase co-existence line 
[21] [25]. To look for any singular behaviour of bulk viscosity ( at the QCD 
phase transition, that might arise due to the contribution from these modes, we 
consider a general Widom -Griffith's form of equation of state, which describes the 
critical behaviour of magnetization in the vicinity of Tc and is compatible with 
the singularities due to these Gold-stone modes. 
It is given by [26j 

y = fix) (27) 
where the variables y and x are given by 

y = h/M\ X = t/M^'^ 
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The variables t and h are the normahzed reduced temperature t = {T — Tc)/Tq 
and the magnetic field h = H/Hq. The normalization of function / are taken to 
be 

/(0) = 1, /(-i)=o 

From the critical exponents (3 and 6 one can derive all other critical exponents. 
The various critical exponents are related as 

di/ = /3(l + <5), 7 = /3(<5-l), iyc = iy/P6 

Another way of expressing the dependence of magnetization M on the reduced 
temperature variable t and magnetic field h is 

M = h^'^fcit/h^l^^) (28) 

where, fa is the scaling function. The scaling forms in Eq.(27) and Eq.(28) are 
clearly equivalent since the variables x and y are related to the scaling function fa 
and its argument by 

y = /G^ ^ = m^nfc'^' 

In presence of Gold-stone modes the scaling function fa is given by [27] 

fa{z) = {-zf{l + ~c,(3i-z)-^'/') (29) 

where z = t/h^/^^ . As has been discussed in [23], Eq.(29) can be easily obtained 
from the magnetic equation of state derived by Wallace and Zia p8] . The universal 
amphtude C2 is given in Table 1. 
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Integrating the expression M ~ dGs/dH, where magnetization M and scahng 
function fo are defined by Eq.(28) and Eq.(29) respectively, the free energy G^, 
comes out to be 



= h{-t)^ + 2/3 C2 /3 i-tf-^'/^ h^''' 



(30) 



To describe the QCD phase transition which is controlled by two independent 
parameters, charge chemical potential // and temperature T, we once again chose to 
parametrize the reduced temperature variable t as, t= t + with t= [T — Tc)/ 
Tc and fi — n/Tc- Using the above form of free energy G^, we will calculate all 
the relevant quantities required for examining the scaling behaviour of C- As in 
previous case we will again limit our study to = case only. 
First we calculate the scaling of susceptibility xtt — d'^Gg/dT'^. Evaluating the 
double derivative of Gg with respect to temperature T , we get 



Using the first and second order derivatives of the reduced temperature variable t, 
dt/dT— T^^ and dH/dT'^— respectively, the scahng form of susceptibihty xtt 
becomes 



Next we calculate the scaling behaviour of susceptibility Xix/j, — 

d^Gs/dfjL^ Defining 
a double derivative of G^ with respect to baryon chemical potential /i, we get 



XTT i-tf-'^ [I + 52{-t)-^'/'] i-dt/dTf 



Xtt ~ {-t) 



13-2 



(31) 



15 



Using the relations dt/dfi = 2Aii/Tc = and dH/d/j,^ = 2A/Tc, the scahng form 
of susceptibihty becomes 

X,, - Hf-' (32) 

Now finally we calculate the scaling behaviour of x^t — d'^Gs/d/idT. Evaluating 
the double derivative of free energy with respect to fj, and T we get. 

X,T ^ i-tf-'^ [l + C2(-t)-^^/2] {dt/dll){dt/dT) 

Using the relations dt/dT = dH/d^dT = and dt/d^ = 2A/xTp^ we get 

X^T = (33) 
Therefore the scaling behaviour of specific heat Cy defined in Eq.(14) is given by 

{-tf-^ (34) 



Cv^T 



xIt 

Xtt 



which for the static critical exponents of 3d 0(2) model, has a divergent behaviour 
near the QCD phase transition. Now as already mentioned in previous section the 
divergent behaviour of bulk viscosity ^ at critical point is controlled by the term 
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Using again Eq.(12) and Eq.(13) together with the scahng laws for susceptibihties 
X/u^i) XnT, Xtt and specific heat Cy we have 

[dP\ _ sx^f, - nXj^T /, N/3-2 



\9nJ^ CvXnn 

Also, using Eq.(28) in conjugation with Eq.(29) we have the scaling relations 
dE/dT = {M/E){dM/dT) ~ anddE/dn = {M/E){dM/dn) ~ {-tf^-\ 

Therefore the singular part of the bulk viscosity ( scales as 

Now using the static critical exponents for 3d 0(2) model, given in Table 1, the 
singular part of the bulk viscosity, C**^, vanishes near the QCD phase transition. 

VI. Results and Discussions. 

We studied the non-equilibrium properties of a quasi-particle medium at finite tem- 
perature and density in a kinetic theory under the relaxation time approximation, 
where the mass of each quasi-particle is assumed to be temperature and chemical 
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potential dependent. In particular we studied the scaling behaviour of bulk viscos- 
ity ( near the QCD chiral phase transition in the 3d 0(2) universality class. It is 
found that the bulk viscosity ( approaches to a finite value, (^''^9"'"'' = (^— ^(^'^^^ = o) , 
near the QCD chiral phase transition. This scaling behaviour of ( prevails even 
after incorporating the effect of Gold-stone modes, that arise due to the explicit 
breaking of the continuous 0(4) symmetry in the staggered fermion formulation of 
the lattice QCD at non-zero lattice spacing. A similar non-divergent result for ( is 
reported in ref [29], where it is found that the the ratio (/s (here s is the entropy 
density) has a maximum value of (/s\max — 0.06 at critical temperature Tc, in 
contradiction to the results reported in ref wherein a sharp rise of fraction (/s 
near the QCD phase transition was obtained, using an analysis based on the exact 
sum rule for QCD with 2 + 1 flavors, lattice QCD data on e — 3P and a sim- 
ple model for the spectral function. This non-divergent behaviour of ( at critical 
temperature Tc using the non-perturbative method of Ads/CFT correspondence 
is similar to leading order perturbative QCD result of (/s ~ 0.02 |3U] for a 
coupling strength of 0.06 < ag < 0.3. On the other hand the Gold-stone modes 
are found to have a significant effect on the scaling behaviour of the specific heat 
'Cy' near the QCD phase transition. These modes lead to the divergent behaviour 
of the specific heat at the critical temperature Tc- However these divergences do 
not appear in the bulk viscosity as the scaling of 'C in our quasi-particle model 
is always determined by the product ~ (M^/Cy) {dM /dT + dM /dfx), conse- 
quently oc Cy^, which is different than the scaling behaviour of 'C given by 
^div ^ (j^^ found in ref [S]. This behaviour of bulk viscosity under the kinetic 
theory in relaxation time approximation should not be regarded as essentially new. 
In ref [31] it is shown that a gas of structure less point particles have negligible 
bulk viscosity both in the non-relativistic and extreme relativistic limits. It is 
particularly due to this reason that mostly bulk viscosity C, is often neglected even 
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when the shear viscosity is taken into account. However it should be mentioned 
that our analysis of scaling behaviour of bulk viscosity ( is valid only under the 
relaxation time approximation. An extension beyond this approximation that al- 
lows for a proper incorporation of the long range fluctuations of the soft modes 
might possibly result in the divergent behaviour of bulk viscosity ( at the QCD 
phase transition, with the divergent behaviour of bulk viscosity being controlled 
by the dynamic critical exponents rather than the static critical exponents used in 
this study. 
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